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We examine the formation and evolution of primordial black holes (PBH's) after hybrid inflation. Our goal is 
to assess the effects of various theoretical uncertainties on the extrapolation from a given inflation model to a 
spectrum of primordial black hole masses. The context of our work is an examination of the possibility^ 
that the dark matter is comprised of Planck- mass black hole remnants (BHR's). As an example we focus on a 
particular scenario|3,l3 m which the black holes form from quantum perturbations that were generated during 
hybrid inflation. We find the correspondence between hybrid inflation parameters and the range of initial PBH 
masses that would allow BHR's to comprise the dark matter, taking account of the possible early presence of 
radiation and its accretion onto the PBH's. 



1. INTRODUCTION 

The main goal of this paper is to assess of the ef- 
fects of various theoretical uncertainties on extrapo- 
lation from a given inflation model to a spectrum of 
primordial black holes. These uncertainties include 
initial size of the PBH's, efficiency of accretion of am- 
bient radiation onto the PBH's, and possible presence 
of radiation due to reheating mechanisms acting prior 
to PH formation. As a particular example we exam- 
ine a scenario 0, H in which: (1) significant numbers 
of PBH's form just after an epoch of hybrid inflation, 
(2) reheating of the Universe occurs primarily through 
Hawking evaporation of these PBH's to Planck scale 
black hole remnants (BHR's), and (3) the dark mat- 
ter consists of these BHR's. In this BHR dark matter 
(BHRDM) scenario the problem is further simplified 
because there is a fairly sharply defined character- 
istic mass scale at which the PBH's form. (BHR's 
as a dark matter candidate were first proposed by 
MacGibbon|a[: other authors have also explored this 
possibility|alg). 

The endstage of Hawking evaporation is not well 
understood. However, heuristic arguments^, Q sug- 
gest that black holes might not evaporate completely, 
but instead may leave behind a stable Planck mass 
remnant, and that a black hole's temperature is given 
by 
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rather than the Hawking value 
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We use Planck units: Planck's constant h, the speed 
of light c, Newton's gravitational constant G, and the 



Boltzmann constant kj, are all set to 1. Thus masses 
are in units of Planck mass mp = y/hc/G ~ 2.2 x 
1CT 5 g (« 1.2 x 10 19 GeV), lengths in units of Planck 
length lp = y/ Gh/c 3 w 1.6 x 10~ 33 cm, times in units 
of Planck time tp = \J Gh/c 5 ~ 5.4 x 10~ 44 sec, etc. 
We use a hybrid inflation^ IIo| potential given by: 

F^^) = (M 2 -^ 2 ) 2 + imV 2 + ^V 2 ■ (3) 

The effective mass of field tp goes from positive to 
negative as <j> decreases from large values. This shift 
to a negative ("tachyonic") mass for ip occurs when 

1/2 



decreases to the critical value 
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There are two inflation regimes in this model: (1) slow 
rolling of <f> down "trough" where <f> > (f> c , with tp w 0, 
(2) Rapid fall of tp to tp±, beginning when (p reaches 
4> c . Large perturbations occuring at the "phase tran- 
sition between these two regimes later produce PBH's 
at a fairly sharply defined mass when they re-enter 
the horizon. 

As discussed in Ref. [TlllT^|. one finds two evolution 
equations for the simple "PBH+radiation epoch" of a 
homogeneous, isotropic, flat, Friedmann Universe: 
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Here a = — , a prime ' means d/da, and we use scaled 
black hole and radiation energy densities defined by: 
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where rij is the initial number density of PBH's, and 
Pbh = riiOT^M is the energy density in PBH's. Note 
that Jbh is also equal to the mass M of the PBH's. 
The first term in the numerator on the right hand side 
of Eq.Q) represents the accretion of radiation (with 
energy density p r ) onto the black holes, and the second 
term repesents Hawking evaporation. Here we've used 
the simpler result in Eq.(J2J), and in simulations dis- 
cussed later will simply put in a cut-off at the Planck 
mass at the end of evaporation. The factor g ~ 100 
gives the multiplicity of particles at the high tempera- 
tures characteristic of black holes near the end of their 
evaporation. In obtaining these equations we used the 
high frequency ("geometrical optics") limit |l3| for the 
absorption cross section of the black holes, which is 
27irM 2 for all relativistic particles. This limit is well- 
satisfied for the situations considered in our scenarios. 
However, for reasons to be discussed later, we have 
also introduced into Eq. [5] an "efficiency factor" for 
accretion, F. 



2. PBH FORMATION 
2.1. "Standard picture" 

We begin by reviewing the treatment of PBH for- 
mation used in Refs. 0, El which is based on the 
Press-Schechter^^ type of formalism put forward by 
Carr|l5j. Assuming spherically symmetric density 
perturbations with Gaussian radial profile and rms 
amplitude 5{M), and an equation of state p = wp with 
Q < w < 1, Carr argued that the probability P(m) of 
a region of mass M collapsing to form a PBH is given 



by 
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Assuming a "hard" equation of state (in our case, p = 
p/3) and a flat Universe the PBH's are expected to 
form with approximately the horizon mass. 

The dependence of initial conditions f B H,i> fr,i ( as 
well as the initial number density of PBH's m appear- 
ing in ct\ and oti) upon the hybrid inflation potential 
of Eq. |3| can be reduced to two parameters and 
s, where H* « y/8iv /3A/ 2 is the Hubble parameter 



during inflation and s = — ~ 



The 



evolution of tp during the waterfall regime is given by: 
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where the subscript "ie" denotes the end of inflation. 

The number of e-foldings of inflation between the 
"phase transition" near <fi ~ <fi c and the end of inflation 
is given by 
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The PBH's form when the fractional density per- 
turbations 5 occurring at the "phase transition" near 
<p ~ <pc reenter the horizon. This density perturbation 
is closely related to the "curvature perturbation" 1Z 
at horizon re-entrv|l6|: 



5 ^Sp = 2 
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Here p — wp is the equation of state, and w at the 
time of re-entry is denoted by w reen t- We also have 



K = 
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where t reen t denotes that the quantity in brackets is 
to be evaluated at the time when the perturbation of 
interest enters the horizon. Furthermore, the average 
spectral perturbation due to quantum fluctuations is 
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If the Universe is approximately matter-dominated 
between the end of inflation and horizon re-entry, then 
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while if it is more nearly radiation-dominated between 
the end of inflation and horizon re-entry, we have 

Hi *\e-™°H* . (14) 



Regardless of where between the two extremes the 
actual evolution falls, we assume that the Universe has 
become radiation dominated by the time of re-entry, 
i.e. w reent = 1/3. 

The standard assumption is that PBH's form with 
the horizon mass. As will be discussed below, this may 
be only roughly true, and changes from this are im- 
portant when considering the effects of accretion onto 
the PBH's. So we introduce a factor xbh(< 1) which 
gives the ratio of the initial PBH Schwarzschild radius 
to the horizon radius H^ 1 at the time of formation: 



f BH ,i = xbhH^/2 



(15) 



Finally, the initial proportions of PBH's and radia- 
tion (where everything except the PBH's falls into the 
category of "radiation") is obtained as follows. From 
Eas. HTirrei and|Slwe have K ~ l/s. Putting this into 
Ea.irni gives 
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Assume the probability P(M) of formation is given by 
Eq.0 where S(M) is given by 8 reent from Eq.fSI (Re- 
call M is the mass inside the horizon at re-entry, and 
thus approximately equal to the mass of the PBH's 
that form.) The initial proportions of PBH's and ra- 
diation is then just 

fr,i _ 1 - P{M)X B H 
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2.2. Uncertainties in PBH formation 
threshold and initial mass 

There remain significant uncertainties in both the 
threshold overdensity for PBH formation and in the 
size of the PBH's that do formal El 111 III EJ 
I22I |23| . Early work by Nadejin, Novikov, and 
Polnarev 20| and Novikov and Polnarev 22] gave re- 
sult that black holes form with order 10% the hori- 
zon mass and also that threshold S is higher. Work 
by some authors |2ll 123] indicates critical phenomena 
in gravitational collapse, leading to the possibility of 
black hole formation with significantly less than the 
horizon mass. Green, et.al.[lj| concluded as the re- 
sults of a study comparing the Press-Schechter based 
formalism with a peaks formalism that the "standard" 
formulation is fairly good if the threshold is in the 
range 0.3 to 0.5. We do not attempt to settle these 
uncertainties here, but note for later use that the ini- 
tial mass of the PBH's may be at least somewhat less 
than the horizon mass. 
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3. SIMULATION AND RESULTS 

Given a choice of the hybrid inflation parameters 
s,-H*, we can calculate initial values / ri ,, /bha and 
then solve Eqns.0]for f r (a) and Jbh (a) up to the time 
when the PBH's have evaporated to remnants. After 
that, f r and Jbh remain constant in time (note in 
particular that Jbh = !)• Matter-radiation equality 
occurs at 
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Here the subscript "eq" denotes matter-radiation 
equality and the subscript "e" denotes the end of PBH 
evaporation. Evolution then continues into the era 
which is dominated first by dark matter and eventu- 
ally, as the present time is approached, by dark energy. 
We assume the dark energy behaves like a cosmologi- 
cal constant, so that fj e = Pde is constant in time. 
(We ignore baryons, in effect just lumping them in 
with the dark matter since they have the same equa- 
tion of state, p w 0.) 

The present age of the universe is given by 



to 



du 
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Since the influence of dark energy is negligible until 
relatively late, the time of matter-radiation equality 
is given by 



teq t e 



du 



+ 



(20) 



Figure 1: Case of no accretion, w — between the end of 
inflation and the time of PBH formation, and xbh = 1 
(baseline case discussed in Refs. 110). The top plot 
shows the line in H, , s parameter space satisfying the 
constraint that matter-radiation equality occurs at the 
observed redshift (z — 3234). (Here we display H, in 
conventional units of GeV rather than Planck units.) 
The next plot covers the same region in the H* , s plane 
and shows the time t e at which the evaporation of the 
PBH's to BHR's is completed; the gray scale coding in 
the latter plot is shown at the bottom of the figure. 



where t e is the PBH evaporation time. 

Provided that matter-radiation equality occurs in 
our simulations at the observed redshift z 3234 (and 
the PBH's do indeed complete their evaporation), the 
constraint that to is the observationally inferred value 
13.5 Gyr is also satisfied. 

3.1. Baseline model 

We begin by showing results for the case discussed 



Refs. 



for which accretion was not included 
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Figure 2: Case with accretion included, w — between 
the end of inflation and the time of PBH formation, 
xbh = 1, and F = 1. The meaning of the plot is as 
described for the top plot in Figure Note that the 
result m this figure cannot be exactly correct, since the 
PBH's initially grow to be larger than the horizon - see 
Figured and discussion in text. 




13 14 15 16 

Log 10 H, 

Figure 4: Loci in H t , s parameter space satisfying the 
constraint that matter-radiation equality occurs at the 
observed redshift (z = 3234), for Xbh ~ 0.4, with 
accretion (solid line) and without accretion (dashed line); 
i.e., same type of plot as in top part of FigureQ Other 
parameters are: w — between the end of inflation and 
the time of PBH formation, and F = 1 for the case with 
accretion. 
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Figure 3: For particular point on the line shown in 
Figure H(s = 2.45, log 10 [H,/GeV] = 14.8), we show the 
PBH radius 2}bh (solid line) and the Hubble horizon 
H^ 1 (dashed line), both normalized by /sj?,i; this is as 
obtained from formulation of accretion used in Eq.[l] 
(with F = 1 in Eq.^J. 



and to which we shall refer as the "baseline case". 
Thus it was assumed that (for whatever reason) the 
"accretion efficiency factor" F is essentially zero. We 
also assume (1) black holes form with horizon size 
{xbh = 1), and (2) w ~ between the end of in- 
flation and the time of PBH formation. The result 
is shown in Figure fl] The top plot shows the line in 
H* , s parameter space satisfying the constraint that 
matter-radiation equality occurs at the observed red- 
shift (z — 3234). The next plot covers the same region 
in the , s plane and shows the time t e at which the 
evaporation of the PBH's to BHR's is completed; the 
gray scale coding in the latter plot is shown at the 
bottom of the figure. A safe criterion would be that 
they have evaporated by the time of the electroweak 
phase transition, expected to be around a TeV, or 
10 -10 sec. The actual values are typically somewhat 
more than this (by one or two orders of magnitude), 
but the criterion is somewhat uncertain anyway. 



3.2. Accretion with F = 1 and x b h = 1 

In Figure|21we include accretion in accordance with 
Eq. 2| (with F = 1 in Eq. |5J|. There is significant 
change from the case of no accretion. Closer exami- 
nation shows that in this formulation, where the black 
holes are assumed to form with the horizon mass (i.e., 
Xbh — 1) an d the accretion efficiency F = 1, the 
black holes accrete essentially all the available radi- 
ation background, before the scale factor grows ap- 
preciably, and the early growth rate of the black holes 
is faster than the growth of the horizon. This is illus- 
trated in Figure |3 which shows (as a function of scale 
factor a) the PBH radius (solid line) and the Hub- 
ble horizon H^ 1 (dashed line) for the case s = 2.45, 
log 10 [H*/GeV] — 14.8 (which is near the middle of 
the line in Figure 

As is obvious from causality (and was long ago 
pointed out by Zeldovich and Novikovpij) the mass 
of the PBH cannot increase faster than the amount of 
mass within the sound horizon after formation, which 
in turn is less than the optical horizon (exactly equal 
to H^ 1 assuming radiation domination, and in any 
case equal up to factor of order 1) 

The faster-than-horizon growth indicates a break- 
down of one or more of the assumptions in our formu- 
lation of accretion, for example, (1) that the radiation 
background remains smooth and uniform, and (2) that 
the black holes form with the horizon mass. 

3.3. Cases with x bh = 0.4 

As was also pointed out by Zeldovich and 
Novikovpij. if PBH's form with significantly less than 
the horizon mass, accretion is unimportant. To get a 
limit on the possible effect of accretion in the BHRDM 
scenario, we reduce Xbh just enough to avoid having 
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Figure 5: The PBH radius 2/bh (solid line) as a function 
of scale factor a, along with the Hubble horizon H^ 1 
(dashed line), both normalized by fBH.i, at s = 3.95, 
log 10 [.ff*/GeV] = 13.85, for x BH = 0.4 (top plot) and 
x bh = 0.41 (bottom plot). 
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Figure 6: Loci in H,,s parameter space satisfying the 
constraint that matter-radiation equality occurs at the 
observed redshift (z = 3234), for accretion with F = 0.4 
(solid line) and without accretion (dashed line, same as 
in Figure 0. 
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PBH's grow to be larger than the Hubble horizon 
H~ l , The required value is Xbh — 0.4, and Figure 01 
shows the loci in H*,s parameter space satisfying the 
constraint that matter-radiation equality occurs at the 
observed redshift (z = 3234), for this case. The result 
with accretion is the solid line, and for comparison 
the dashed line shows the result without accretion. 
The effect of accretion is greatest at larger s (upper 
part of plots) because the parameter y = f r ,i/ fsH.i is 
larger, i.e. there is initially more radiation available 
for accretion. In Figure [S] we illustrate the growth of 
the PBH's (solid line) and the growth of the horizon 
(dashed line) for a case near the top of the solid line 
in FigureH at s = 3.95, log l0 [#*/GeV] = 13.85. The 
top plot is for Xbh — 0.4; here the PBH's gain signif- 
icant mass by accretion, but never grow faster than 
the horizon. The bottom plot is for Xbh = 0.41; here 
the PBH's do exceed H^ 1 for a short time. 

3.4. Reduction of accretion efficiency 

The solid line in Figure^] represents only an upper 
limit on the PBH growth for this x = 0.4 case, as there 
is no guarantee that our accretion formulation is valid 
even though the growth is not faster than the horizon. 
For example, it could well be the case that the as- 
sumption that the radiation background remains uni- 
form is still violated in reality, i.e. the region around 




Figure 7: The PBH radius 2}bh (solid line) as a function 
of scale factor a, along with the Hubble horizon H -1 
(dashed line), both normalized by }bh,i, at s = 3.95, 
\og 10 [H t /GeV] = 13.85, for F = 0.4 (top plot) and 
F = 0.41 (bottom plot). 



the black holes could become depleted of radiation, 
making F < 1. 

It is of course possible to avoid the faster-than- 
horizon growth in the Xbh — 1 case as well, if the 
accretion efficiency is reduced. Reduction to F = 0.4 
is sufficient, and results are shown in Figures and 
(We would expect F — ► 1 eventually, but by that time 
accretion becomes neglibible anyway.) 

We note Hacyan|23 found (in a general-relativistic 
treatment based on Einstein- Strauss vacuole model) 
that the initial growth of PBH's by radiation accretion 
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Figure 8: Locus in _ff*,s parameter space satisfying the 
constraint that matter-radiation equality occurs at the 
observed redshift (z = 3234), for case with no accretion 
included, w — 1/3 between the end of inflation and the 
time of PBH formation, and xbh = 1- 



could be near the horizon rate if the initial PBH size 
near horizon size, but also interpreted this result as a 
generous upper limit on growth rate. 

3.5. Possible additional reheating 
mechanisms 

The baseline case assumed that the Universe is ap- 
proximately matter dominated from the end of infla- 
tion until nearly the formation time of the PBH's. 
However PBH evaporation is not necessarily the only 
reheating mechanism. We need not exclude the possi- 
bility of additional reheating mechanism(s) that could 
result in approximate radiation domination well be- 
fore PBH formation. Thus we compare also a case 
(Figure |SJ where w = 1/3 between the end of infla- 
tion and the PBH formation time. Here we assume 
other parameters are as in Figure F°r a given ini- 
tial s, -ff*, the w = 1/3 case gives smaller initial mass 
of PBH's than w = 0, since M oc 2e 2N " instead of 
M cx |e 3A,c . 



4. SUMMARY AND CONCLUSIONS 

We examined effects associated with early presence 
of radiation upon extrapolation from hybrid inflation 
parameters to a spectrum of PBH's in the BHRDM 
model. It appears that accretion can produce some 
early rapid growth. This may occur even if the re- 
quirement that the PBH's not grow to be larger than 
the horizon H~ x at any time is enforced by reducing 
either the fraction Xbh of the horizon mass going into 
the initial PBH mass or by reducing the accretion ef- 
ficiency F. However our results should be taken as an 
upper limit on the effects of accretion, as there is no 
guarantee that the efficiency factor is close to F = 1 
during the early rapid-accretion phase. Both these ac- 
cretion effects and possible radiation production well 



before PBH formation affect the correspondence be- 
tween the hybrid inflation parameters s, and initial 
PBH masses, but do not seriously change the results 
of the BHRDM scenario. 
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